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Abstract 

We consider the relativistic Vlasov-Maxwell and Vlasov-Nordstrom systems which describe 
large particle ensembles interacting by either electromagnetic fields or a relativistic scalar 
gravity model. For both systems we derive a radiation formula analogous to the Einstein 
quadrupole formula in general relativity. 
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1 Introduction and Main Results 

This paper is an investigation of the mathematical properties of certain models for the interaction 
of matter, described by a kinetic equation, with radiation, described by hyperbolic equations. The 
first model, the relativistic Vlasov-Maxwell system, plays an important role in plasma physics. The 
motivation for studying the second model, the Vlasov-Nordstrom system, comes from the theory 
of gravitation. On a mathematical level the Vlasov-Maxwell system can also give insights into 
gravity. 

The most precise existing theory of gravitation, general relativity, predicts that certain as- 
trophysical systems, such as colliding black holes or neutron stars, will give rise to gravitational 
radiation. There is a major international effort under way to detect these gravitational waves jH]. 
In order to relate the general theory to predictions of what the detectors will see it is necessary 
to use approximation methods - the exact theory is too complicated. The mathematical status of 
these approximations remains unclear although partial results exist. This paper is intended as a 
contribution to understanding the mathematical structures involved. 

* Supported in parts by DFG priority research program SPP 1095 



Since the solutions of the equations of general relativity are so difficult to analyze rigorously it is 
useful to start with model problems. One possibility is the scalar theory of gravitation considered 
here, the Vlasov-Nordstrom theory [Hj. It has already been used as a model problem for numerical 
relativity in [2*Tj . 

Among the approximation methods used to study gravitational radiation those which are most 
accessible mathematically are the post-Newtonian approximations. Some information on these has 
been obtained in jTTj and [T%j . Results which are analogous to these but go much further have 
been obtained for the Vlasov-Maxwell and Vlasov-Nordstrom systems in [3] and [2] respectively. 
None of these results include radiation explicitly. Here we take a first step in doing so. On the 
other hand, for the case of finite particle systems interacting with their self-induced fields there 
are several rigorous results concerning radiation; see (22] for an up-to-date review. 

Our main results (Theorem 11.41 and Theorem 11.91 below) are relations between the motion of 
matter and the radiation flux at infinity for the Vlasov-Maxwell and Vlasov-Nordstrom systems 
respectively. They are analogues of the Einstein quadrupole formula |2*3*1 (4.5.13)] which is a basic 
tool in computing the flux of gravitational waves from a given source. In the case of the Einstein 
and Maxwell equations a spherically symmetric system does not radiate. For the Vlasov-Nordstrom 
system a spherical system can radiate and the specialization of the general formula to that case is 
computed. In [23] a difference between the spherically symmetric and the general case was claimed 
but we have not succeeded in connecting this to our results. The main theorems are obtained under 
plausible assumptions on the behavior of global solutions of the relevant system (Assumption 11.11 
and Assumption II .61 below). The former can be proved to hold in the case of small data. 

For the systems we are going to consider the (scalar) energy density e and the (vector) momen- 
tum density V are related by the conservation law 

d t e + \7-V = 0. 

Defining the local energy in the ball of radius r > as 



£r(t) = / e(t,x)dx, 

J\x\<r 

this conservation law and the divergence theorem imply that 
d 



S r (t)= d t e(t,x)dx = - V-V(t,x)dx = - x-V(t,x)da(x), (1.1) 

^ J \x\<r J \x\<.r J \x\=r 

where x — o denotes the outer unit normal. More specifically, for the relativistic Vlasov-Maxwell 
system with two particle species, 

e R yu(t,x)=c 2 J ^/l + c-^(f + + f-)(t,x,p)dp+^-(\E(t,x)\ 2 + \B(t,x)\ 2 ), (1.2) 

V RVM (t,x) =c 2 J p(f + + f-)(t,x,p)dp+^-E(t,x)xB(t,x), (1.3) 

whereas for the Vlasov-Nordstrom system, 

e YN (t,x)=c 2 J v /l + c -V/(t,a;,p)dp+^((^(t,x)) 2 + c 2 |V0(t,x)| 2 ), (1.4) 

,4 



f C 

Vy N (t,x) =c 2 pf(t,x,p)dp- — d t (f)(t,x)V<f)(t,x 



2 



Our assumptions on the support of the distribution function will be such that the contributions of 
/ p(f + + f )dp to Prvm and / pf dp to Pvn vanish for \x\ —r large. Hence we arrive at 

l^ RVM (t) = ^_ f x-(BxE)(t,x)da(x) 
at An J\ x \ =r 

for the relativistic Vlasov- Maxwell system, and 

^-£^{t) = ^- / x-(dt<l>V<t>)(t,x)d<T(x) 
dt 4irJ\ x \ =r 

for the Vlasov-Nordstrom system. 

The main results of this paper are concerned with the expansion of these energy fluxes for 
r,c — s>oo and \t — c _1 r| < const. Under suitable assumptions we will prove that, to leading order, 

±s^(t)^-±\dMu)\ 2 , 

where u = t — c~ 1 r denotes the retarded time and T>(u) = J xpo(u,x)dx is the dipole moment asso- 
ciated to the Newtonian limit of the relativistic Vlasov- Maxwell system. Similarly, 

with a more complicated radiation term 1Z associated to the Newtonian limit of the Vlasov- 
Nordstrom system. In the spherically symmetric case, d t 7Z(u),u) is found to be proportional 
to d t £^i n (u), the change of kinetic energy of the Newtonian system. The exact statements are 
contained in Theorems 11.41 and 11.91 below. 



1.1 Dipole Radiation in the Relativistic Vlasov-Maxwell System 



The relativistic Vlasov-Maxwell system describes a large ensemble of particles which move at pos- 
sibly relativistic speeds and interact only by the electromagnetic fields which the ensemble creates 
collectively. Collisions among the particles are assumed to be sufficiently rare to be neglected [T3] . 
In order to see effects due to radiation damping it is necessary that there are at least two species 
of particles with different charge-to-mass ratios. For the sake of simplicity we assume that there 
are exactly two species with their masses normalized to unity and their charges normalized to 
plus and minus unity, respectively. The density of the positively and negatively charged particles 
in phase space is given by the non-negative distribution functions f ± = f ± (t,x,p), depending on 
time tGM, position a;GlR 3 , and momentum pGM 3 . Their dynamics is governed by the relativistic 
Vlasov-Maxwell system 



9t/ ± +p ■ V s / ± ± (E + c~ x p x B) ■ V p / ± = 0, 
cVxE = -d t B, cVxB = d t E+4iTj, 
V-E = 4np, V-5 = 0, 

p= [ (f + -r)d P , 3= f p(f + -r)d P , 



(RVMc) 



where 



P = 1P, 7 ; 



;i + c-Y)- l/a , p 2 = H 2 , and 1=1 

J JR3 



;i.5) 
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The electric field E = E(t,x) G M 3 and the magnetic field B = B(t,x) G M 3 satisfy the wave equations 



-d 2 + c 2 A)E = 4n(c 2 Vp + d t j) and {-d 2 +c 2 A)B = -AixcV x j. 



;i.6) 



In order to determine the radiation of the system at infinity, we have to consider solutions that 
are isolated from incoming radiation. For the wave equations in this means that we need to 

restrict ourselves to the retarded part of the solutions. Accordingly, (|RVMc|) is replaced by 



dt^ +P- V*/* ±(E + c-'p x B) ■ V p / ± -- 
E(t,x) = -J ( Vp + c- 2 d tJ ) (t - c- 1 \y - x\ ,y) 



= 0, 

dy 
y-x\ 



B(t,x) = c- 1 [ \7xj(t-c- l \y-xly)-^- 

J \y-x\ 

p=j(f + -r)dp, j=jp(f + -f-)dp, 



(retRVMc) 



which we call the retarded relativistic Vlasov-Maxwell system. We prescribe initial data 

f ± (0,x,p) = f ± >°(x,p), x, P eR 3 , 



;i.7) 



for the densities at t = 0; these data do not depend on c. However, the corresponding solution 
(f + ,f~ ,E,B) does depend on c, but we do not make explicit this dependence through our notation. 
We refer to Remark ll.5f c) below for the case of initial data varying with c. Our standing assumption 
is that the initial data are non-negative, smooth, and compactly supported, 



°° m3 xlR 3 ) 



/ ± '° G C 

and we fix positive constants R ,P ,S such that 



r' o >o, 



±,o 



/ '°(x,p) = for \x\>R or \p\>Po, and ||/ 
Every solution of ()retRVMc|) satisfies the identity 

f ± (t,x,p) = f ± >°(X ± (0,t,x,p),P ± (0,t,x,p)), 
where s h- > (X ± (s,t,x,p),P :iz ( y s,t,x,p)) solves the characteristic system 

x = p, p = ±(E + c~ x px B), 



\w 3 -> 



;i.9) 



;i.io) 



with data X ± {t 1 t,x 1 p) =x and P ± {t 1 t 1 x 1 p) —p. Hence < f ± {t 1 x 1 p) < W^^W^- In order to derive 
our results on radiation, we have to assume certain a priori bounds on the corresponding solutions 
of (jretRVMc)) . In particular, the latter have to exist globally in time. 



Assumption 1.1 (a) For each c>l the system AretRVAlb]) has a unique solution f^EC 2 



x 



), EeC 



2/D xl 3 : 



), BeC 



2/D xl 3 : 



, satisfying the initial condition \1. 



(b) There exists Pi>0 such that f ± (t,x,p) = for \p\> P\ and all c>l. In particular, / ± (t,a;,p) 
for \x\ > Rq + Pi \t\ by HTTP . 
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(c) For every T>0, R>0, and P>0 there exists a constant Mi(T,R,P)>0 such that 



\d? +1 f ± (t,x,p)\ + \d?VJ ± (t,x,p)\<M 1 (T,R,P) 
for \t\<T , |x|<P ; |p|<P, and a = 0,1, uniformly in c>l. 
Note that none of the constants in Assumption 11.11 may depend on c. The constants 

Ro>Po,So,Pi,Mi 

from (jl.fjj) and Assumption 11.11 are considered to be the "basic" ones. Any other constant which 
appears in an estimate is only allowed to depend on these. Checking the arguments from [3 IB], it 
can be shown that Assumption II. II holds at least for sufficiently "small" initial data / ±,c . A more 
precise investigation of the set of initial data leading to solutions which satisfy Assumption 11.11 is 
not part of this paper. The main point we want to make here is that whenever Assumption 11.11 is 
verified, then the technique described below can be employed. 

We will need estimates relating the solutions of (jretRVMc)) to the corresponding Newtonian 
problem obtained in the limit c — >oo. This sort of information usually goes under the name of 
post-Newtonian approximation; see jT3 E] . For this, one formally expands the solutions in powers 
of c -1 as 

/ ± = /o ± + c- 1 A ± + c- 2 / 2 ± + -, 
E = E + C - 1 E 1 + C - 2 E 2 + ..., 
B = B + c- 1 B 1 + C - 2 B 2 + ..., 

with coefficient functions f^ 1 , Ej, and Bj independent of c. Moreover, by (|1.5j) . 

P = P- (c- 2 /2)p 2 p + . . . , 7 = 1- (c- 2 /2)p 2 + .... 

These expansions can be substituted into (jretRVMc|) . and comparing coefficients at every order 
gives a sequence of equations for the coefficients. The Newtonian limit of (jretRVMc|) is given by 
the plasma physics case of the Vlasov-Poisson system: 

d t f±+p-V x f±±E -V p f± = 0, 



E (t,x) 



x-y 



p {t,y)dy, 



\x-y\ 

(fo-fo) d Pi 
ft(0,x,p) = f ± >°(x,p) 



(VPpl) 



The following proposition addresses the well-known solvability properties of ( VPpl ). Clearly, 
(fo,fo,E ) is independent of c, and we refer to e.g. [2011111 for the regularity of the solution. 



Proposition 1.2 There are constants R2,P2>0, and for every T>0, R>0, and P>0, there is a 
constant M 2 (T,R,P) >0, with the following properties. For initial data f ± '° as above, there exists 
a unique global solution (f^,E ) of \VPpHjj so that 

(a) /^^(Ixtfxl 3 ) and P e x M 3 ;M 3 ) ; 
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(b) if \t\<l, then f^(t,x,p) = for \x\>R 2 or\p\>P 2 , 

(c) if\t\<T, \x\<R, \p\<P, and a = 0,1, then 

| d? fo(t,x,p)\ + \d?E (t,x)\<M 2 (T, R,P). 

For the approximation of solutions of (jretRVMc)) by solutions of QVPpip , we state the following 
result without proof; the result follows like the analogous one for (jRVMc|) . cf. [T§1 H]. 

Proposition 1.3 Choose the constants P\>0 and Mi(T,R,P)>0 according to Assumption II. 11 
Then for every T>0, R>0, and P>0 there are constants M 3 (T,R,P)>0 and M±(T,R)>0 
with the following property. If c>2P 1; let (f ± ,E,B) and (/ ,E ) denote the global solutions of 
HretRVMc]) and \VPpHjj provided by Assumvtion M . 1\ and Proposition U .2l respectively, with initial 
data as above. Then 

(a) \f ± (t,x,p)-f^(t,x,p)\<M 3 (T,R,P)c- 2 for\t\<T, \x\<R, and \p\<P, 

(b) \E{t,x) - E {t,x)\<M 4 (T,R)c- 2 for \t\<T and \x\<R, 

(c) \B(t,x) \ <M 4 (T,R)c~ 1 for \t\<T and \x\<R. 

It is important to note that all the "derived" constants R 2 ,P 2 ,M 2 ,M3,M4 appearing above do only 
depend on the basic constants R ,P ,S ,Pi,Mi. We are now ready to state our first main result. 

Theorem 1.4 (Radiation for (jretRVMcp ^ Put r* = max{2( J R + P 1 ),R 2 } and 

M RVM = {(t,r,c):r>2r*,c>2P 1 , |t-c -1 r| < 1, r > c 3 }. 

If (t, r, c) G A^RVM? then with r=\x\, x= and u = t — c~ 1 \x\, 

x-(Bx E)(t,x) + c~ 4 r- 2 \x x d 2 V{u)\ 2 < A(c' 5 r~ 2 + c~ 2 r~ 3 + c~ 1 r~ 4 ), (1.12) 
for a constant A>0 depending only on Ro,P ,S ,Pi,Mi. In particular, 

d Sf YM {t) = ^-\ x-(BxE)(t,x)da(x) 



dt r w 4tt 



|a;|=r 



= -^\d 2 V(u)\ 2 + 0{c^ + c- 1 r- 1 +r- 2 ) (1.13) 
6c 6 

for (t,r,c)eM K vM- Here £™ M (t) = f lx ^ r e RVM (t,x)dx, see (OJ), and 

V{u)= J xp (u,x)dx 
denotes the dipole moment associated to the Vlasov-Poisson system \VPpHjj . 

Remark 1.5 (a) The condition r>c 3 in TWrvm is not needed for the proof of ()1.12j) and ()1.13|) . 

It just guarantees that c~ 2 r~ 3 <c~ 5 r~ 2 and c _1 r _1 <c -4 . 

(b) The same estimate ()1.12j) can be derived, possibly with a different constant A, if the condition 
\u\ < 1 is replaced by \u\ <uq for some constant uq > 0. 
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(c) As long as the constants Rq,P ,Sq,Pi,Mi remain independent of c, one can also allow for c- 
dependent initial data ff' , both for (jretRVMc)) and QVPpip . However, in this case the functions 
(Jo ,Eq) become c-dependent, too. For instance, in the particular case 

/ c ± ' = /o ± '° + c- 1 /i ± '° + c- 2 /5 , 

with /5 t ' ,/{ t '°, and f^° satisfying suitable bounds (independently of c for f^°), Theorem 11.41 
remains valid, if L and E are replaced by the approximations / + c~ 1 / 1 and Eq + c^Ei, re- 
spectively. Here (f^,E ) is the solution of QVPpip for the initial data fo'°, and (f^,Ei) solves the 
Vlasov-Poisson system linearized about (f^,E ), under the initial condition /^(O) = 

(d) In the case of one species only, say / _,o = 0, there is no dipole radiation, since then <9 2 D = 0, 
cf. (l2~T3j) below. 

(e) For spherically symmetric solutions there is again no dipole radiation. In fact, if po(t,—x) = 
Po(t,x) for 16K 3 , then T> = by symmetry. 



The proof of Theorem 11.41 is given in Section 12.11 



1.2 Monopole Radiation in the Vlasov-Nordstrom System 

If we set all physical constants (except the speed of light c) equal to unity, then the Vlasov- 
Nordstrom system is given by 



dtf+p-VJ- [OS0)p + c 2 7 V0j • V P / = 4(S0)/, 
(-<9 t 2 + c 2 A)0 = 47r/i, 



ifdp, 



(VNc) 



where we continue to use the notation from p.5jl . and where S = dt+p-'V. The matter distribution 
is modeled through the nonnegative density function f = f(t,x,p), whereas the scalar function 
4> = <f)(t,x) describes the gravitational field. We refer to [HI El H E] for the global existence of 
smooth solutions to (jVNcjl . In analogy to the passage from (jRVMc|) to (jretRVMcj) . the solutions 
of (|VNc|) that are isolated from incoming radiation are the solutions of the retarded system 



dtf+p-VJ- (50)p + c 2 7 V0 ■V p f = 4(S<f>)f, 



<j)(t,x) 



p(t-c \y-x\,y)- 



dy 



\y-x\ 



ifdp, 



(ret VNc) 



which we call the retarded Vlasov-Nordstrom system. We continue to make the standing hypotheses 
JEBJ) and (EHl) for the initial data f(0,x,p)=f°(x,p) of (jretVNcjl . A solution of IJretVNcjl satisfies 
the relation 

f(t,x,p) = f o (X(0,t,x, P ),P(0,t,x,p))e^\ (1.14) 
where si— > (X(s,t, x,p),P(s,t,x,p)) denotes the solution of the characteristic system 



x—p, p— —(Scf))p — c 2 'y'V(f), 



1.15) 
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with X(t,t,x,p) — x and P(t,t,x,p) —p. This implies that as long as the solution exists 



o<f(t,x }P )<\\r\Li 

note that 0<O. Concerning solutions of (jret VNcjl . we make the following 

Assumption 1.6 (a) For each c> 1 the system hetVNc]) has a unique solution f EC 2 (RxM 3 x 



>eC 2 



, satisfying the initial condition f(0,x,p) = f°(x,p). 



(b) There exists Pi>0 such that f(t,x,p) = for \p\>Pi and all c>l; by M.15)) this 
implies that f(t,x,p) = for \x\ ^-Ro + 'Pil^l- 

(c) For every T>0, R>0, and P>0 there exists a constant Mi(T,R,P) >0 such that 

for \t\ <T, \x\ <R, \p\ <P , and a= 1,2. In addition, for every T>0 and R>0 there exists 
a constant Mi(T,R) >0 such that 

\ ( t>{t,x)\ + \V(t ) {t,x)\ + \d t <P{t,x)\<M l {T,R) 

for \t\ < T and \x\ < R, uniformly in c> 1. 

Again Rq,P ,So,Pi,Mi are considered to be the "basic" constants, all other constants being derived 
from these. We remark that for "small" initial data the existence of global-in-time solutions is 
shown in jl2j, where also bounds on the solutions are obtained. It is reasonable to expect that these 
solutions have the required regularity for smooth initial data, cf. [10], and that on compact time 
intervals estimates as in Assumption II .61 (c) can be derived uniformly in c. The crucial assumption 
is the bound on the momentum support in part (b), which needs to be uniform in c as well. 

The Newtonian approximation for c — > oo of (jretVNc|) is found by means of the formal expansion 

f = fo + c- 1 f 1 + c~ 2 f 2 + ..., 

= 00 + C~ Vl + C~ 2 02 + C" 3 3 + C" 4 04 + • • • , 

see [TO] 12] • Thereby it is verified that this (lowest order) Newtonian approximation of (|retVNc| ) is 
given by the gravitational case of the Vlasov-Poisson system 

$/o+P-V b /o-V02-V p /o = O, 

Po(t,y) 



Po 



x-y\ 
fodp, 



■dy, 



(VPgr) 



f {0,x,p) = f°{x,p). 

The analogue of Proposition 11.21 is valid for (JVPgrJ). Note that (fo,4>2) is independent of c. 

Proposition 1.7 There are constants R2,P2>0, and for every T>0, R>0, and P>0, there is 
a constant M2(T,R,P) >0, with the following properties. For initial data f° as above, there exists 
a unique global solution (/o,02) of \VPgr\j so that 
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(a) f eC°°(RxR 3 xR 3 ) and <p 2 € C°°{R x R 3 ), 

(b) if\t\<l, then f (t,x,p) = for |x|>-R 2 or\p\>P 2 , 

(c) if\t\<T, \x\<R, \p\<P, and a = 0,1,2, then 

\d?fo(t,x,p)\ + \dt +1 Ut^)\ + \d?VMt,x)\ < M 2 (T,R,P). 
By Pj, we also have the following rigorous result concerning the Newtonian limit of (jretVNcjl . 

Proposition 1.8 Choose the constants P\>0 and Mi(T,R,P)>0 according to Assumption \1.(A 
Then for every T>0, R>0, and P>0 there are constants M 3 (T,R,P)>0 and M 4 (T,i?)>0 with 
the following properties. If c>2P±, let (f,(j)) and (f 0,^)2) denote the global solutions of \retVNc\l 
and \VPgr\) provided by Assumption \l.b\ and Proposition respectively, with initial data as 
above. Then 

(a) \f(t,x,p)-f (t,x,p)\<M 3 (T J R,P)c- 2 for\t\<T, \x\<R, and \p\<P, 

(b) |V0(t,x)|<M 4 (T, J R)c- 2 for \t\<T and \x\<R, 

(c) \d t <p(t,x)-c- 2 d t <p 2 (t,x)\ + \V<p(t,x)-c- 2 VM t ^)\< M 4( T ,R)c~ i for \t\<T and \x\<R. 
After these preparations we can state our second main result. 

Theorem 1.9 (Radiation for (jretVJNc^ Put r* = max{2(i2 + P 1 ),i? 2 } and 

Mv N = {(t,r,c):r>2r*,c>2P 1 , |t — c _1 r| <l,r>c 6 }. 
If (t,r,c) G .MvN; then with r= \x\, ^= jfj; and u = t — c~ 1 \x\, 

x-(d t (j)V(f))(t,x) + c^ 9 r- 2 (d t n(x,u)) 2 <A(c- 10 r- 2 + C - 4 r- 3 ), (1.16) 
for a constant A>0 depending only on Rq,Pq,Pi,Mi,Sq. In particular, 

d £™(f) = / x-{d t (f>V(f>){t,x)da{x) 



dt r y ' 4vr 



|x|=r 



l — [ {d t TZ(to,u)) 2 da(to) + 0(c- & + r- 1 ) (1.17) 



47TC 5 

for (t,r,c) G M-vjsi- Here £^ N (t) = f,, <r ey^(t,x)dx, see \l.4\l , and 

TZ(x,u) = -^ J \x-V(j) 2 {u,y)\ 2 dy- J j (x-p) 2 f (u,y,p)dpdy + 4:£ kin ( y u), (1.18) 

where 

S kin (t) = ^ J J p 2 f Q (t,x,p)dpdx (1.19) 
denotes the kinetic energy associated to the Vlasov-Poisson system \ VPgr\) . 

Defining £ po t(t) = — ^ / 1 V(fi 2 (t,x)\ 2 dx, the total energy £ (t) = £ kin (t) + £ pot (t) is conserved along 
solutions of dVPgr |. 



Remark 1.10 (a) Once again the condition r>c 6 in .Mvn is not needed for the proof of 

It only has to be included in order that the second error term 0(c~ A r~ 3 ) is at least as good as the 

first one, which is O(c~ 10 r~ 2 ). 

(b) In the sense of Remark |l .51 (b) and (c), one could allow for \u\ <u and/or c-dependent initial 
data. 

For spherically symmetric solutions, Theorem 11.91 simplifies as follows. 

Corollary 1.11 (Radiation for spherically symmetric solutions to (JretVNcj)) 

Define = max{2(R Q + Pi),R 2 } and 



M 



VN : 



j(t,r,c) : r>2r*,c>2P 1 , \t- c -1 r| < 1, r > c 6 |. 



If (t,r,c) G .Mvn; then with r = \x\ and u = t — c V , 

(«9 t <R0) (t,x) + ^ c" 9 r" 2 {d t S Un (u)) 2 1 < A(c- W r- 2 + c^r- 3 ) , 
for a constant A>0 depending only on Rq,Pq,So,Pi,Mi. In particular, 

±£™(t) = ^J _ {d t <t>d r <t>){t,x)da{x) = -^{d t 8^ n {u)) +0(c- 6 + r 
for (t,r,c)eMvj<i- 

The proofs of Theorem 11.91 and Corollary 11.111 are carried out in Section 12.21 

2 Proofs 

2.1 Proof of Theorem 11.41 

To expand E(t,x) and B(t,x) as given by (|retRVMc|) . we recall from Assumption 11.11 (b) that 
f ± (t,x,p) = for |x|>i2o + Pi|t|. It follows that p(t,x) = and j(t,x) = for \x\ >R + Pi\t\. If 
(t,x,c) G Mrvm, then \u\ = \ t — c _1 |x| | < 1 and c> 2P X . Thus if \y \ >2(R + Pi), then 

-Ro + -Pi|^ — c_1 |l/ — x \\ — Ro + -Pi|m + c _1 |x| — c _1 |?/ — x|| 

< J Ro + Pi(kl+c~ 1 |y|)<i?o + Pi(l + (2Pi)- 1 |y|)<|y|. 

Hence _F(t — c _1 |?/ — = for both F — —(Vp + c~ 2 dtj) or F = c _1 V x j. Thus for the ^/-integrals 
defining E and S in (jretRVMc)) . it is sufficient to extend these over the ball \y\ <max{2(i? + 
Pi),R 2 } = r*. 

In what follows g = 0(c~ k r~ l ) denotes a function such that 

\g(t,x)\ < Ac~ k r~ l for all \x\=r>2r*, c>2P 1 , and \t — c _1 |x| | < 1, 

with A only depending on the basic constants. The following lemma states a representation for E 
and B similar to the Friedlander radiation field; see [T5J p. 91/92] and jHj- 
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Lemma 2.1 The fields can be written as 

E{t,x) = E r&A (t,x) + 0{r- 2 ) and B{t,x) = B™ d {t,x) + 0{ C - 1 r- 2 ), 

where 

E™ d (t,x) = -r' 1 [ {Vp + c- 2 d t j){u + c- l x-y,y)dy, (2.1) 
J\y\<r* 

B rad {t,x) = c^r- 1 Vxjiu + c^x-y^dy. (2.2) 



Proof: Consider E first, and let F = — (Vp + c 2 d t j). According to Assumption II . II (c) . we have 
< AMi(l+r*,r*,p*) — 0(1) for some constant A>0, where = max {Pi,P 2 } and 



(,..) = (t-c 1 \y-x\,y) = (u + c 1 \x\-c 1 \y-x\,y). 
N 

i ui _ u. _ 

r _1 + 0(r- 2 ) 



If Ixl = r > 2r* and lyl < r*, then r^-r < H < 2. It follows that 

ii — * ici i — \v—x\ — \x — r* — 



1 




1 


- + - 


x\ 


-\y 


— x 




\y-x\ 




x\ 


\y-x\ 


\x\ 





for all \y \ <r*. Therefore by (jretRVMc)) . 

E(t,x)= [F(...) T ^L_= [ F(...) r ^-r= / FG.Op + ^r- 2 ))^ 
= r- x [ F(...)dy + 0{r~ 2 ). 

J\y\<r, 

Next we note that for \y\ <r* and \x\ —r >2r*, 



\x 



— \x — y\ = \x\ — \x\ y/l — 2x- y/\x\ + \y\ 2 /\x 

1 



r\ -\.,-\ | l + -(-2x-y/\x\ + \y\ 2 /\x\ 2 )+0{r- 2 ] 
= x-y + 0(r-- 1 ). (2.3) 

Since 

\F(...)-F(u + c- 1 x-y,y)\<\\d t F\\ L00 c- 1 \\x\-\y-x\-x-y\ = 0(c- 1 r- 1 ) 

by Assumption 11.11 (c) and (|2.3|) . we get i? = £' rad + (9(r _2 ). The proof for the magnetic field is 
analogous, using F = c^ 1 V X j. □ 

Now we need to investigate the relation between E ra,d and B raA . For this, we recall the continuity 
equation d t p + V ■ j = and calculate 

Vp(*) = V, [/>(*)] +c~ 1 x V, ■ [?'(*)] - c~ 2 (x ■ dtj(*))x, 
Vx j(*) = V y x [j(*)]-c" 1 xx<9 4 j(*), 

where 

(*) = (u + c _1 x -y,y) 
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is the argument. This follows just from evaluating the total derivatives. Since L y \ <r dy = J dy in 
(12. 1|) and ()2.2|) by Assumption 11.11 (b), integration by parts shows that all V^-terms drop out. 
Consequently due to u = t — c~ 1 r the relations 

E™ d {t,x) = -r- 1 [ [Vp(*)+c- 2 d tJ (*)}dy 
J\v\<r* 

•~' ' -c~ 2 (x-d t j(*))x + c~ 2 d t j(*) dy 

jiu + c^x-y^-ix-jiu + c^x-y^^x 



\y\<r* 



-c 2 r 1 d i 



B rad (t,x) = c~V 



\y\<r* 

Vxj(*)dy 



dy, (2.4) 



y\<r* 



= — c 2 r 1 d t / xxj(u+c l x-y,y)dy 
J \y\<r* 

are obtained. Note that in particular £' rad and B Tad are of the same order in c _1 and r _1 , i.e., 

E vad (t,x)=B vad (t,x) = 0(c- 2 r- 1 ) (2.5) 
by Assumption 11.11 (c). Observing 



x x 



\y\<r* 



j(*)-{x-j(*))x 



dy= I xxj(*)dy, 
\y\<r* 



differentiation w.r. to t yields the important formula 

xx£ rad (^) = B rad (t,4 

Also 



x- 



y\<r* 



j(*)-(x-j(*))x 



dy = 0, 



so that 



x-E™ d (t,x)=x-B™ d (t,x) = 0. 
Collecting the results from Lemma [2. II and (J2.5|) . it follows that 

x ■ (B x E) = x ■ ([5 rad + C(c-V- 2 )] x [E Tad + 0{r- 2 )} 

= x ■ (V ad x £ rad + 0{c~ 2 r- 3 ) + 0(c- 3 r- 3 ) + Oic^r^ 
= -\xx E rad \ 2 + 0(c- 2 r- 3 ) + C(c- 1 r- 4 ), 

since by flUHJ) and (|2~7I) . 
Eqns. (Q and (J23J) imply 



x • (5 rad x £ rad ) =i-([k £ rad ] x £ rad ) = x ■ ({x ■ £ rad ) E iad - | £ rad 1 2 x) 

rad 1 2 I X X E rS,d 1 2 



x ■ (B x E) = —c r 



4„-2 



ix (9 t / j('« + c %-y,y)dy +0(c r 3 ) + 0(c r 

\y\<r* 



(2.6) 



(2.7) 



(2- 



(2.9) 
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To expand the square as c— >oo, we note that \p\ >p*>Pi implies / ± (t,y,p) = for all £eR and 
all t/Gl 3 by Assumption 11.11 (b). Therefore we can always replace the average over momentum 
space J dp by J lpl < p ^dp. For \p\ <p*, 

V p p = 7id M 3 - c~ 2 rfp®p = id R3 + 0(c~ 2 ) 
by p.5jl . Furthermore, using Assumption II. 1| 

v a ,(/ + -r)(*) = v»[(/ + -r)(*)]-c- 1 «ft(/ + -/-)(*) 

= V tf [(/ + -/-)(*)]+l { | y |< r ., H < P .}0(c- 1 ). 

Utilizing this, (jretRVMcJ) . and Proposition IO[ we get, writing (*,p) = (u + c -1 a; ■?/,?/, p), 

dtj(*)dy= / pd t (f + -f~)(*,p)dpdy 
y\<r* J J 



P^-P-v x (f + -r)(*,p) 

-{E + c~ l p x 5) • V p (/ + + D(*,p)^ dp<fy 
0(c~ l ) + ! ( V p p(E + c- 1 pxB)(f + + f-)(*,p)dpdy 

J\y\<r* J\p\<p* 

0(c- 1 ) 

+ / / (id R s + O(c- 2 ))(E + O(c- 2 ))(f+ + f +O(c- 2 ))(*,p)dpdy 
J\y\<r*J\v\<P* 

I I E (f+ + U)(*,p)dpdy + O(c- 1 ). (2.10) 



Also 



by Proposition O(c). Thus and (jlHTTjl yield 



2 



x-( J BxE) = -c- 4 r- 2 xx / / E (f+ + fo)(u,y,p)dpdy+O( C - 1 ) 

l^l^: r « J \p\^zP* 

+ 0{c- 2 r- 3 ) + 0{c' 1 r- 4 ) 

= -c~ 4 r~ 2 xx J J E (f£ + fo)(u,y,p)dpdy 

+0(c~ 5 r' 2 ) + 0(c- 2 r- 3 ) + 0(c' 1 r- 4 ), (2.11) 

since by Proposition 11.21 (b). /q (u,y 7 p) — for \y\>r*> R 2 or \p\>p*> Pi- Defining the dipole 
moment 

with po from ( |VPplD , we obtain by the Vlasov equation in ( |VPplD that 

d t V = j J x d t (f+-f )dpdx = - J jx(p-V x (f+-U) + E -V p (f+ + U))dpd. 



X 



J j P(fo~fo )dpdx 
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and 

d 2 t V = j j p d t (f+-f,)dpdx = - j fp(p-V x (f+-fv) + E .V p (f++fo))dpdx 

= j j E (f+ + f )dpdx. (2.12) 

Due to (j2~TT[) it follows that 

x-(BxE) = -c- 4 r- 2 \x x d 2 V{u) \ 2 + 0(c~ 5 r - 2 ) + C(c^ 2 r" 3 ) + C^c" 1 ^ 4 ), 
which completes the proof of ([1.120 . Concerning ([1.13(1 . we have 

d S Rv M{t)= c f x-(BxE)(t,x)da(x)-c 2 [ [ (x -p){f + + f-)(t,x,p)dpda(x) 

dt *XJ\x\=r J\x\=rJ 

by p. 1(1 and ([1.3(1 . For \x\ — r and (t,r,c) G-Mrvm, we can estimate Rq + P\\t\ = R + P\\u + c~ l r \ < 
R + P\ + \ < r = \x\, since r >2r*>2(R + P 1 ). Therefore f ± (t,x,p) = by Assumption Q (b) , 
and this yields 



dt 



47T J|x|=r 



|x|=r 

Hence for (j!.13() it suffices to use ([1.12(1 and to note that 

-3^-2 

l\x\=r ^f^J\u,\=l 

by integration. □ 



■^-J \xxd 2 t V{u)\ 2 da{x) = -^j Juxd 2 V(u)\ 2 da(uj) = -^\d 2 V(u)\ 2 



Proof of Remark ESJ (d): If / ~(t = 0) = /" ° = 0, then also / ~ = by (0711). Thus defining 
(j) (t,x)=ff\x-y l ^fo (t,y,p)dpdy = J\x-y\ ~ l p (t, y) dy, we get £ O = -V0 O and A0o = -4vrp o . 
Consequently by ([2.12(1 . 

d f 2 £> = J j E f+dpdx = J E p dx= j ' V<f) o A(f) o dx = 0. (2.13) 

Hence there is no dipole radiation in this case. □ 

2.2 Proof of Theorem fCT 

By (|retVNc|) . d t (f> and V</> are given by 



d t (f){t,x) = -c 2 J dtn{t-c x \y-x\,y) V 
V(j)(t,x) = -c~ 2 J Vii{t-c- l \y-x\,y)-^- 
In full analogy to Lemma 12. 1( we obtain the following representation. 



\y-x\ 
dy_ 

x\ 
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Lemma 2.2 We can write 

d t (f)(t,x) = (<9 t 0) rad (t,x) + 0{c' 2 r- 2 ) and V<p{t,x) = (V0) rad (t,x) + 0( C - 2 r - 2 ), 

where 

\rad/ + „\ „-2„-l / a ../„. , -1 



(dtcjyy^x) = -c- 2 r- L / d tf M(u+c- l x-y,y)dy, (2.14) 

■/|j/|<r» 

(V0) md (t,x) = -c _2 r _1 / V^{u + c' l x-y,y)dy. (2.15) 

■/ |2/]<r* 

Let again (*) = (m + c~ 1 x •?/,?/) denote the argument. Then V y [/i(*)] =c~ 1 xd t fi(*) + V/x(*). Since 
Jj^l <r cfo/ = J in (|2.15j) . it follows that 

x ■ ( V0) rad (t,x) = -c^r' 1 x ■ J V//(*) dy = c^r^x ■ J xd t /i{*) dy = -c _1 (<9 t 0) rad (t, x) . 

The same argument shows that (V0) rad = 0(c~ 3 r -1 ), and also (<9 t 0) rad = C(c~ 2 r _1 ) by ()2.14j) . 
Hence we find from Lemma [2,21 that 



x-(d t (j)'V(j))(t,x) = —c r 



5^,-2 



d t fi{*)dy 



2 

4^-3 



+ C( C - 4 r- 3 ). (2.16) 



In order to expand the square we use, following [M], the differential operators 

T = c~ 1 xd t + V and S = d t + p-V. 

Then 

9 t = (l-c" 1 p-x)" 1 (5-p-T) 

and V y [/!(*)] = T/i(*) is a total derivative. Hence the corresponding term drops out upon integra- 
tion with respect to y. Observing the relation 

V p -[(50)p + c 2 7 V0]=3(50), 

the Vlasov equation in (jretVNc|) yields 

)dy = J J 'jd t f(*,p)dpdy 

= J J^ii-c-'p-xy'^s^p+c^v^-vj+Ais^fy^^dpdy 



V p (7(l-c- 1 p-x)^ 1 )-[(50)p + c 2 7 V0]/(*,p)^ o ?|/ 
+ J J 1 {l-c- 1 p-xY\S ( t>)f{* 1 p)dpdy, (2.17) 

where i*,p) = {u + c~ 1 x-y 1 y,p). A direct calculation shows that 

V p (7 (1 - c^p-xY 1 ) = V p (( v 7 ! +P 2 /c 2 - c _1 p - =7 2 (l-c- 1 p-x)~ 2 ( C - 1 S-c- 2 p). 
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If \p\ >p*>Pi, then f{*,p) = by Assumption 11.61 (b). Furthermore, if \u\ <1 and \x\ >2r*, then 
\y\ >r* enforces /(*,£>) = as before. Therefore we can replace JJ dpdy by f\ y i <r $\ p \< p dydp in 
the integrals occurring in ()2.17|) . In other words, we may always assume that both \y\ and \p\ are 
bounded, with a bound depending only on the basic constants. Accordingly, 

7 = 1 + C(c- 2 ), 7 2 = 1 + C( C - 2 ), p = 1P = p + 0(c- 2 ), (2.18) 
(l-c- 1 p-x)- 1 = l + 0(c- 1 ), (l-c- 1 p-x)- 2 = l + 2c- 1 p-x + 0(c- 2 ). (2.19) 

This results in 

V p (7 (1 - c^p ■ x)- 1 ) = c~ l x + cT 2 (2(p • x)x -p) + 0(cT 3 ). (2.20) 
Furthermore, since |w + c -1 x -?/| < l + r#, also 

f(*,p) = f (*,p) + O(c- 2 ), 
(S<P)(*,p) = c~ 2 (S ( p 2 )(*,p) + 0(c- i ) = 0(c- 2 ), 
V0(*) = c - 2 V0 2 (*) + O(c- 4 ), 

by Proposition 11.81 and (|2.18|) . where S4>2 = d t 4>2+p-'V4>2- Observe that here the constants M 3 (l + 
r*,r # ,p*) and M4(l+r*,r#) enter the bounds on C(c~ 2 ) and (9(c -4 ). Hence from (|2.17|) . (|2.18j) . 
(EUnj), and dZZBl we get 



dtfi{*)dy = - / (c _1 x + c _2 (2(p-x)x-p) + C?(c- 3 )) 

•%!<»•» J\p\<p, 

■ O(c- 2 ) + (l + O(c- 2 ))(v0 2 + O(c- 2 ))] (/o(*,p) + 0(c- 2 ))^ 
+ /,|<, /,<, ( 1 + ( C " 2 ))( 1 + °( C " 1 ))( C_2 (^) + C, ( C " 4 )) 

7o(*,p)+0(c- 2 )Wi/ 



-C X' 



{l + 2c- l p-x)V{(f)2fo)(*,p)dpdy 

y\<r*J\p\<p* 



+c~ 2 [ [ (Sfo+P-V^foi^dpdy + Oic- 3 ). (2.21) 

./|j/|<r» ^|p|<p* 

Let ^ denote either V02 or <9 4 02- Then by Proposition 11.71 (c). 

(V'/oX*,?) = (V ; /o)(M + C~ 1 X-t/,?/,p) 

= (V'/o) («,2/,p) + c -1 (x • y) (V»/ ) (u,y,p) + 0{c~ 2 ) 
= (ijf )(u,y,p) + O(c- 1 ). 
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Hence fl2~2T|) yields 



d t ^*)dy = -c- 1 x- I / {(V(l> 2 fo)(u,y,p)+c- 1 (x-y)d t (V ( f )2 f )(u,y,p)+0(c- 2 ) 

\y\<r* J |p|<p* 



+2c 1 (x-p)V(f) 2 f (u,y,p))dpdy 
+c- 2 I I (Sfo+ P -V4>2)fo(u,y,p)dpdy + O(c- 3 ) 

J \v\<r* J p|<r>* 



y\<r* J\p\<p* 

(x-V<p 2 )po(u,y)dy 



J k 



\y\<r* 

-2t ' 



-c d t / (x-y)(x'Vfa)po(u,y)dy 

J\v\<r* 

+c~ 2 [ [ (S(P 2 + p-V<j )2 -2(x-p)(x-V ( j )2 ))fo(u,y,p)dpdy } (2.22) 
J\y\<r* J\v\<p* 

recalling that u = t — c~ 1 r. In view of Proposition 11.71 (b) we may extend all integrals over the 
whole space again. Now 

V(j) 2 p {u,y)dy = J j V(f) 2 (u,y)f (u,y,p)dpdy = 

by Lemma 1231 below, whence the lowest order term drops out. In addition, Lemma 1231 also shows 
that 

r r _ 

(S(j) 2 )(u,y)fo(u,y,p)dpdy = -2dtS kin (u) 



as well as 

/• / 

(p-V(j) 2 )(u,y)fo(u,y,p)dpdy = -d t £ kin (u) 



and 

(x-p)(x-V(f) 2 ){u,y)fo(u,y,p)dpdy = ~d t J j (x-p) 2 f (u,y,p)dpdy. 
Finally, we can also write 

(x-y)(x-V(f) 2 )(u,y)po{u,y)dy = -£p t{u)--^ J \x-V(j) 2 {u,y)\ 2 dy 

by Lemma l2~3l and since \x\ = 1. Using this and d t £ pat = —d t £km (see the remarks following (|1.19p ) 



in (|2.22J) . and collecting all the terms, it follows that 

d t ^(*)dy=-c- 2 d t n(x,u) + 0(c- 3 ), (2.23) 

with 1Z(x,u) as in (jl.l8|) . Inserting (|2.23|) into ()2.16|) . we see that 

x-(d t (j)W(j))(t,x) = -c- 5 r- 2 |-c- 2 a t ^(^,u) + C(c- 3 )| 2 + C(c- 4 r~ 3 ) 
= -c- 9 r- 2 (d t TZ(x,u)) 2 + O(c~ 10 r- 2 ) + C( C - 4 r" 3 ). 

Therefore 16|) is proved. Concerning (jl.l7j) . the fact that 



d .£^{t) = ^- I X'(O l oY<>)(t..r)<l*(.r) 

' \x\=r 
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dt r y ' 4vr 



is due to (t,r,c) G A^vn, analogously to the argument in the proof of Theorem II .41 Hence (|1.17|) is 
a direct consequence of (jl.l6j) . changing variables as x = ru. □ 

We still need to give the proof of 
Lemma 2.3 For the Vlasov-Poisson system \VPgr\ j, 

(S(f) 2 )(t,x) f (t,x,p) dpdx = -2d t £km(t) , 
J J p-V(p 2 (t,x)fo(t,x,p)dpdx = -dt£km(t), 

J J ^•p)^-V(h)(t,x)f (t,x, P )dpdx=~d t j J (i-pfh(t,x,p)dpdx (£gM 3 ), 

J j ^■x)^-V<p 2 )(t,x)f (t,x,p)dpdx = -\^\ 2 £ pot (t)-^ J |e-V0 2 (t,x)| 2 rfx (£eM 3 ), 

where £ kin (t) = ^ ff p 2 f (t,x,p)dpdx, see M.lty . and £ pot (t) = —-^ f \V(j) 2 (t,x)\ 2 dx. 
Proof: Firstly, since A02 = 47rpo, 

V<f) 2 (t,x)f (t,x,p)dpdx = / V<f) 2 (t,x)p (t,x)dx = — jV<j) 2 (t,x) A<p 2 (t,x)dx 



^ / V-\VMt,x)\ 2 dx = 0. 
on J 



For the remaining assertions we define the mass current density as jo — / pfodp. Integration of the 
Vlasov equation with respect to p implies the continuity equation d t pa + V ■ Jo — 0- Hence 



S(f) 2 f dpdx = J (d t (f) 2 p + V(f) 2 -jo)dx = J (dtfapo-fo V- j )dx 
= d tf (p2podx = 2d t £pot(t) = -2d t £ kin (t) 



by conservation of energy, and 



p-Vfafodpdx = I j -V(f) 2 dx= I d t p (j) 2 dx 

1 



\x-y\ 
d t £ pot (t) = -d t £ kiQ (t). 



-d t po(t,x)p (t,y)dxdy 



Furthermore, by flVPgr ) , 



dtjj (^p) 2 f dpdx = J j (e-p) 2 V p -(V0 2 / o )^ 
= ~ 2 J J(Z-P)(t-V<h)fodpdx. 
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For the last assertion, using A0 2 = 47rp o , 

J V<p 2 )p dx = -|- J {i-x)i i d i (j)' 2 d j d j 4)2dx 

= -g^/ (^-V|V0 2 | 2 d:r-^ y |£- V0 2 | 2 dx 

= |^/ |V0 2 | 2 rfx-^ y |£-V0 2 | 2 cfe 

= -|^ 2 ^pot-^ y K-V0 2 | 2 dx. 
This completes the proof of the lemma. □ 

2.3 Proof of Corollary flTlTl 

In this section we verify Corollary 11.111 by specializing Theorem 11.91 to spherically symmetric 
functions. We recall that initial data f° are said to be spherically symmetric, if 

f (Ax,Ap) = f (x,p) 

for any matrix A G S0(3). Then the solution (fo,<p 2 ) of QVPgrD provided by Proposition II . 71 remains 
spherically symmetric for all times. Therefore 

f (t,Ax,Ap) = f (t,x,p), p (t,Ax)=p (t,x), and fa(t,x) = <fr 2 (t, Ax) (2.24) 

holds for all AG SO (3). 

Firstly, this implies V <p 2 = xd T (f) 2 as well as |V0 2 | 2 = \d r <j) 2 \ 2 , d r denoting the radial derivative. 
By choosing AgS0(3) such that Ax = ej (the j's unit vector in M 3 ), (|2.24j) yields 



x- V<p 2 {u,y)\ 2 dy = -j- / \dj<f) 2 {u,y)\ 2 dy = -j- I j^d r <f) 2 {u,y) 



An An. Air 



\y\ 



2 

dy 



If 2 

\d r 4> 2 (u,y)\ 2 dy = --£ pot (u). 



UttJ 1 v " 3 
Similarly, (|2.24|) and |x| 2 = l implies that 



y y (x-p) 2 f (u,y,p)dpdy = J j p) f (u,y,p)dpdy=^ J j p 2 f (u,y,p)dpdy 



= -£kin(«). 

Therefore by ()1.18|) . 

H{x,u) = —^ J \x-V(j) 2 (u,y)\ 2 dy- J j (x-p) 2 f (u,y,p)dpdy + A£ khl (u) 

2 n t \ 10 n / \ 
= -t pot {U) + —t^ m [u). 

Since d t £ pot — — dtS^ by conservation of energy, we get d t TZ(x,u) — ^dtS^^u). Hence Corol- 
lary HHU follows from (ITTnT) and (fTTTTjl . ' □ 
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